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In this article, we calculate the dressed quark propagator with the flat bottom potential in 
the framework of the rain-bow Schwinger-Dyson equation. Then based on the nonperturbative 
dressed quark propagator, we calculate the ir decay constant and the quark condensate. The ir 
decay constant is an important parameter in describing the interplay between dynamical symmetry 
breaking and confinement, while the quark condensate is an order parameter for dynamical chiral 
symmetry breaking. To implement confinement, we prove that the dressed quark propagator has no 
poles on the real timelike p 2 axial, the absence of Kallen-Lehmann spectral representation obviously 
■ precludes the existence of free quarks. 
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I. INTRODUCTION 

The standard model gives highly successful descriptions of hadronic physics in terms of the basic force. Yet its 
key parameters, the masses of the quarks and leptons are specified by the coupling of the Higgs bosons, which is 
I , undetermined. It is possible that the dynamics of the fundamental gauge theories themselves generate the masses 
of all the matter fields, or in other word, dynamical chiral symmetry breaking (for quarks), which corresponding 
^vq , to the absence of the low mass scalar partners of the ir jij. Quantum chromodynamics (QCD) is the appropriate 
theory for describing the strong interaction at high energy. However, the lack of observation of free quarks and gluons 
| leads to the assumption of colour confinement that only quantities which transform as colour singlets can be physical 
observable. One may expect that the nonperturbative infrared behaviour of gluon propagator which may lead to a 

Ph! 

i linear or oscillator potential implements confinement, however, till now, all theoretical studies based on QCD can not 
give satisfactory explanation. Hence we take an equivalent statement that the propagator of a colored state should 
. £h ' have no singularities on the real time-like p 2 axial, the absence of Kallen-Lehmann spectral representation precludes 
' the existence of free quarks 

In fact, dynamical symmetry breaking and confinement are two crucial features of QCD, however, they correspond 

_ _ i 

to two very different energy scales. The possible interplay of the two dynamics lies in the following two facts, the 
first one is that on large energy scale the chiral massive quark interacts with antiquark via a confining potential and 
creates a new quark mass, constituent quark mass; while the second one can be attributed to the double role of the 
7r, as both Nambu-Goldstone boson and qq bound state. 

The Schwinger-Dyson equation, in effect the functional Euler-Lagrange equation of quantum field theory, provides 
a natural framework for investigating the nonperturbative properties of quark and gluon Green's functions. By 
studying the evolution behaviour and analytic structure of the dressed quark propagator, one can obtain valuable 
information about the dynamical symmetry breaking phenomenon and confinement. 
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Although lattice calculations are rigorous in view of QCD, they suffer from lattice artifacts and uncertainties, 
such as Gribov copies, boundary conditions and so on. Furthermore, current technique can not give reliable result 
below lGeV, where the most interesting and novel behaviour is expected to lie. 

It would seem therefore that at the present time Schwinger-Dyson (SD) equation is the most reliable tool for 
studying the infrared behaviour of the quarks and gluons propagators in the continuum limit, while it has its own 
shortcomings. 

The flat-bottom potential (FBP) is a sum of Yukawa potentials, which not only satisfies gauge invariance, chiral 
invariancc and fully rclativistic covariance, but also suppresses the singular point which the Yukawa potential has. It 
works well in understanding the meson structure, such as electromagnetic form factor, radius, decay constant Q. 

Therefore, it is interesting to examine dynamical chiral symmetry breaking and quark confinement with the FBP. 

The article is arranged as follows: in section 1, introduction; in section 2, infrared behaviour of the gluon 
propagator; in section 3, dynamical chiral symmetry breaking; in section 4, quark confinement; in section 5, flat 
bottom potential and Schwinger-Dyson equation; in section 6, conclusion and discussion. 



II. INFRARED BEHAVIOUR OF GLUON PROPAGATOR 



The infrared structure of the gluon propagator has important implication for quark confinement. One might 
expect that the behaviour of the quark interaction in region of small spacelike p 2 determines the long range properties 
of the qq potential and hence implements confinement. 

From the discussion of previous section, we can see that the SD equation is the ideal tool for studying the infrared 
behaviour of the gluon propagator. Here we make a brief skeleton of the results based on Landau gauge and axial 
gauge studies. 

It is impossible to solve an infinity series of coupled SD equations for gluons and quarks (ghosts), so the authors 
have to make truncations both in Landau gauge and axial gauge. In those studies, the quark loop contributions 
are neglected and higher order vertexes are taken as bare or dressed based on the analysis of corresponding Ward- 
Takahashi Identity. 

In axial gauge studies, one makes further assumption that the dressed gluon propagator has the same tensor 
structure as the free one. That is in the following equations let F% — ||, 

D^ip, 7) = Fi(p, 7 )M^(p, n) + F 2 (p, j)N^ u (p, n), 

Mftvtj), n) = dpy - - + n — , 

p ■ n (p ■ ny 

Nftvip, n) = 6^ 5-, (1) 

here 7 is gauge parameter and n is an arbitrary four vector. After solving gluon SD equation one get the conclusion 
that the gluon propagator may have a double pole at the origin, implying an area law behaviour of the Wilson 
loop, which is often regarded as a signal of confinement Further more , the Fourier transform of the time-time 
component of the gluon propagator is related to the one-gluon exchange contribution to the potential between static 
colour charge and a -K, infrared behaviour does generate a linearly rising potential at large distance Q. However, a 
general studies of the properties of the spectral representation of the gluon propagator in axial gauge show that the 
coefficient of the 5^ term cannot be more singular then || . This suggests that neglecting the F2 term is a poor 
approximation because there maybe cancellation in the infrared between F\ and F 2 - 

In Landau gauge studies . Suppose neglecting ghost loop will not remove the dominant infrared contributions 
in QCD also leads to a ^ behaviour for the gluon propagator in the origin. Other studies suggest that fermion loop 
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may suppress the infrared singularity in the gluon propagator, the double pole found in the pure-gauge sector may be 
removed by the nonperturbative fermion loop corrections |lfj,niln7[ . Recently, the gluon SD equation studies including 



ghost loop result in a infrared vanished gluon propagator and a strong infrared enhanced ghost propagator |g4 25 
while lattice calculations support this type conclusion ptjfl . In fact, the truncations applied in these studies may 
impair their persuasion. Furthermore, in some studies of this type, unphysical particle-like singularities occur in the 
dressed ghost-gluon and 3-gluon vertices and the infrared vanished gluon propagator can neither produce dynamical 
chiral symmetry breaking nor compliment confinement p4| , p7| . 

We can take a simple example by analyzing the renormalization group equation for the coupling constant. There 
exists a competing mechanism between the self-interaction of gauge boson and the interaction of boson-fermion, which 
accounts for the different evolution behaviours of the coupling constants with momentum for QCD and QED [jj"2|] . So 
neglecting one another may lead to very different results. 

All theoretical studies based on QCD are suffered from shortcomings in one or another way |Tl| , so they can not 
give satisfactory infrared behaviour for gluon propagator which may lead to a linear or oscillator potential implements 
confinement. 

Here we take an alternate description of confinement based on the analytic structure of the dressed quark prop- 
agator, if the quark Green's function is free of singularities on the real timelike p 2 axial , then denies the existence of 
a single particle spectrum and precludes free quarks. Although by no means a sufficient condition for the formulation 
of finite size qq states, it is an essential ingredient in describing of confinement based on a bold assumption that 
interaction between dressed quarks which does not rise to infinity at large distance. The interaction between dressed 
quarks may diminish at separation beyond the characteristic length scale at which point the quark Green's function 
approaches its vacuum values, but the vacuum's repulsive force due to the absence of a mass pole in the dressed quark 
propagator again compresses the quarks together. 

So the traditional model of confinement based on an effective rising potential can be interpreted as the result 
of combination of the repulsive interaction between quark-vacuum and the attraction of dressed quark-quark with 
increasing separation. Although there maybe a long way to go before the repulsive quark-vacuum interaction can be 
quantified, this description has proven successful in the formulation of a confining nontopological soliton model for 
baryons [Q. In this article, the quark propagator calculated with the FBP does free singularities on the real timelike 
p 2 axial and implements quark confinement. 



III. DYNAMICAL CHIRAL SYMMETRY BREAKING 



The ultraviolet stability of QCD at large spacelike momentum makes a straightforward definition of the chiral limit 
possible, in a manifest contrast to strong coupling QED, whose rigorous definition remains an instructive challenge. 
The 7r decay constant f n is an important parameter in describing the interplay between dynamical symmetry breaking 
and confinement, while the quark condensate serves as an order parameter for dynamical chiral symmetry breaking. 
In the chiral limit, the it plays a double role both as Nambu-Goldstone boson and qq bound state, which lead to the 
result that the scalar part of the inverse dressed quark propagator coincides with the approximated 7r's Bethe-Salpeter 
function. In the following, we calculate the n decay constant , the quark condensate and examine the dynamical chiral 
symmetry breaking phenomenon. 
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A. 7r Decay Constant 



Neglecting ghost, the Ward-Takahashi Identity (WTI) for quark-gluon vertex in the presence explicit quark 



masses can be written as 23 



MWp) = - IS-'ip'hs +1,S-\ P )] + 2imT 5 (p' lP ). (2) 



In the limit of exact chiral symmetry, the WTI reduces to 

k^(p',p) k -^-T l5 B{p 2 ). (3) 

Note that the usual perturbative isovector axial vector vertex is §7^,75. 

In the chiral limit and take k — > 0, the axial vector quark vertex becomes completely dominated by the pseu- 
doscalar coupling of a massless ir to the quark and the subsequent weak decay of the ir into an axial vector current. 
In fact, the mass shell residue of the axial vector quark vertex (T^) is the ir Bethe-Salpeter (r w ) amplitude, in the 
chiral limit can provide a direct connection to the quark propagator JTJ]. The function T^(p+ ■§,£>— §) (r 5 ) contains 
many Dirac structures, such as 757^ (75), 757/i7 • k (757 • k), 757 • pp^ (757 • pp ■ k), ■ ■ • and regular term, singular 
term. Exclude the singular term (T„) , the other terms are supposed to give small contributions and neglected. If 
we solve the Schwinger-Dyson equation for the vertex, we can see their relative weights. In the chiral limit, the 75 
component of r„. makes a good approximation, though it has many Dirac structures, such as 75, 757 • k and so on. 
The 7T decay constant / w is defined by the axial transition amplitude for the on-shell 7r, 

(0\A^{0)\ir n (k)) =iUk fi 5 mn , (4) 

When k — > 0, we obtain 



^ iT^(j>',p)^i{0\A^ m (-k)) ^° -r 5 (p',p)f^. (5) 



Multiplying Eq.(5) with fc M on both sides, 

k ll r%( P ',p) = iUT 5 (p / , P ). (6) 

When combined with Eq.(3), we obtain the Goldberger-Treiman relation for the quark-pseudoscalar vertex, 

r , , . B(pi) . A{p 2 )m{p 2 ) 

r 5 (p, P) = «T 75 — — = iT 75 . (7) 

The simplest generalization of Eq.(7) is 

. B(p' 2 ) + B(p 2 ) 

r 5 (p,p) = "-75 — . (8) 

Based on above discussion, it is possible to obtain an expression for f 2 from the integral equation for the ir decay 
constant. 

z<0|A™(0)k"(fc)> =-j ^tr{iS(p + k)iT?(p + k,p)iS(p)i^-w 5 }. (9) 

As usual, the factor —1 on the right hand side arises from the quark loop. We need to take the axial- vector vertex 
as a perturbative one in order to avoid double counting. The expression for f 2 can be written as 

, f d A p B(p 2 ) ra/ , x , 2v p 2 , dA(p 2 ) dB (p 2 ) . , , . 

f 2 = -12i / , TT, 9 L, o^ {A{p )B{p 2 ) + ^ B-V - A ,V )>- 10 
7 (2tt) 4 (A(p 2 )p 2 - B 2 (jp 2 )) 2 1 v ; v ; 2 V dp 2 dp 2 n K ' 
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In the limit A = 1, we can recover the Pagels and Stokar formula fllSfl . Our numerical results based on the quark 
propagator obtained from the SD equation studies with the FBP do give satisfactory result = 94MeV, we can 
thus fit the parameters in the FBP. 

There is another expression for the n decay constant which determined from global colour model [E6U29] . 



o ,.00 

ft = ^J dssB ( s f{°l - 2 M + sa v a' v ] - s[a s a': a' 2 } s 2 [o"X' - a' v 2 ]} 



(11) 



where 



A(s) B(s) 
v sA(s) 2 +B(s) 2 ' s sA(s) 2 + B(s) 2 ' { ' 

In our calculation, the value of A(p) is far from 1, we prefer Eq.(lO) rather than Eq.(ll), as Eq.(ll) is exact only in 
the limit A(p) = l . 



B. Quark Condensate 

The quark propagator is defined as 

S(x) = (0\T[q(x)q(0)]\0)- (13) 

where q(x) is the quark field and T the time-ordering operator. For the physical vacuum consists of both perturbative 
and nonperturbative parts, so the quark propagator S(x) can be divided into a perturbative and a nonpcrturbative 
part as the following: 

S(x)=Spt{x)+S np (x). (14) 

In the nonperturbative vacuum, the normal-ordered product Snp(x) does not vanish. The widely used nonlocal quark 
condensate (0| : q(x)q(0) : |0) is given by the scalar part of the Fourier transformed inverse quark propagator, 

(0| : q(x)q(0) : 10), = (-4N C ) [ (15) 

In the limit x — > 0, we can obtain the expression for the local quark condensate, 

(01 : mm : 10), = [ dss ^^ . (16) 

In Ref. |l!| the nonlocal condensate is put in the following form : 

<0| : q(x)q(0) : |0) - g(x 2 )(0\ : q(0)q(0) : |0), (17) 

where g(x 2 ) is the vacuum non- locality of the nonlocal quark condensate. The nonlocal properties of the vacuum 
condensates are of principal importance in the study of the distributions of quarks and gluons in hadrons [ p2[ . 
Physically, it means that the vacuum condensates of quark and gluon can flow through the vacuum with nonzero 
momentum. 

The effects from hard gluonic radiative corrections to the quark propagator are connected to a possible change 
of the renormalization scale \i at which the condensates are defined. Those effects are of minor importance for the 
nonperturbative effects in the low and medium energy regions, and neglected in our studies . In this article, we 
take the rain-bow SD equation, so it is not a renormalizable interaction. As isolated closed quark loop integral, the 
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quark condensate has only one energy scale, the momentum cut off, /x, at which the quark condensate is defined. In 
our studies, the energy scale /z is implicitly determined by the effective gluon propagator, we can use Eq.(16) as the 
definition for the quark condensate. However, presently, the nonperturbative technique can not prove the relation 

P " dsS sA^+BH s ) ~ ( L °9(» 2 /A 2 QC n)f, (18) 
at low energy scale. Here d is the anomalous dimension and takes the value d = „„ „ — . 

OJ 33 — 2n-f 

If we take the limit fj, — > oo and assume ultraviolet dominating, the above relation (qq) ^ ~ (Log(/i 2 / AQ CD )) d is 
indeed the case, although the quark condensate is only defined at low energy scale. 

Operator product expansion has proven that at large Euclidean momentum, the effective quark mass evolves as 

Login 2 /A 2 QC D ) 



M-Q 2 h^oo = -^{Log{Q 2 /k 2 QCD )} d - 1 + ™(M){ L J ( y /A ry , (19) 



here c = —^ttj — f\ Vk-^ — vR i s som e constant independent of a |59SJ|, this implies that (01 : qq : |0}m ~ 
[Logi/i 2 1 Aq C D )] d , but not the definition of Eq.(16) . 

On the other hand, the evolution of quark condensate with fj- can be determined from Gell-Mann-Oakes-Renner 
relation and renormalization group equation, 

If we take into account the hard radiative corrections to the quark propagator , there is a more nature definition for 
the quark condensate (in Euclidean space time) |^8| : 

: qq : 0) M = {- — ^ -} -— / d p (21) 

Log{A z JA z QCD ) 16tt 4 J A(p J )V+-B(p ) 

In the second definition, we can take the cut off be A — > cxd, however, at large momentum, the quark propagator 
approximates its asymptotic free form which should be subtracted. From Fig.l, we can see that the cut off A = 2GeV 
makes a sound approximation. 

In this article, fi is taken to be 1 GeV. If we ignore the hard radiative corrections to the quark propagator, from 
Eq.(16), the numerical result gives the quark condensate (0 |gg |0) i = — 218MeV. If we take into account the hard 
radiative corrections, from Eq.(21), we obtain: 

A Q cd = 150MeV, (g(0)g(0))s = -255Me7; (22) 
Aqcd = 200MeF, (q(0)q(0))i = -253MeV; 
K Q cd = 250MeF, (q(0)q(0))i = -25lMeV. 

All our numerical results are compatible with the results of the QCD sum rule approach, and large enough to give 
dynamical chiral symmetry breaking. The definition of Eq.(21) is better. 

IV. QUARK CONFINEMENT 

To demonstrate that quark confinement arises from the absence of singularities on the real timelike p 2 axial, let 
us first make an observation on the large Euclidean time behaviour of the free fermion propagator pcfl . 

f Pr^'^r^i (23) 

v ; J (2ir) 4 k 2 + m 2 y ' 
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The spatially averaged Schwinger function is a particularly insightful tool. Consider the fermion propagator and 
let T = 1x4 represent Euclidean time, then 

<r s (T) = J d 3 x^tr D S(x, T) = \er mT . (24) 

Hence the free fermion's mass can easily be obtained from the large T behaviour of the spatial average: 

lim logcr s (T) = -mT. (25) 

T — >oo 

This is just the approach used to determine the bound state masses in simulations of Lattice QCD. 

In order to demonstrate that the confinement of quarks, we have to study the mass function of the quark and 
prove that there no poles on the real timelike like p 2 axial. So it is necessary to perform an analytic continuation of the 
dressed quark propagator from Euclidean space into Minkowski space p^ — > ipo. However, we have no knowledge of 
the singularity structure of quark propagator in the whole complex plane. One can take an alternative safe procedure 
and stay completely in Euclidean space avoiding analytic continuations of the dressed propagators p2| . Again we 
take the Fourier transform with respect to the Euclidean time T. It is sufficient to consider the scalar part S s , 



S* S (T)= J +C0 ^±e^ T S s 



50 m 

, 2?r q 2 A 2 {q 2 ) + B 2 {p 2 )' { ' 

If S(p) had a pole at p 2 = —m 2 , the Fourier transformed S*(T) would fall off as e~ mT for large T or logS 1 * = —mT. 

In our calculation, for large T, the values of S* is negative, except occasionally a very very small fraction positive 
values which can be safely neglected. We can express S* as |5*|e mir , n is an odd integer, logS* = log \S*\ + inn = 
—mT. If we neglect the imaginary part, we find that when the Euclidean time T is large, there indeed exists a crudely 
approximated linear function with respect to T, which is shown in Fig. 2. However, such a behaviour is hard to acquire 
physical explanation. Here the word 'crudely' should be understand in the linearly fitted sense, to be exact, there is 
not linear function. This precludes the existence of free quarks. 

V. FLAT BOTTOM POTENTIAL AND SCHWINGER- DYSON EQUATION 

The FBP is a sum of Yukawa potentials which is an analogous to the exchange of a series of particles and ghosts 
with different masses (Euclidean Form), 

n 

where N stands for the minimum value of the mass, p is their mass difference, and a,j is their relative coupling 
constant. In Fig. 3, we can see the curve of G{p 2 ) is a gauss-like distribution. Due to the particular condition we take 
for the FBP, there is no divergence in solving the SD equation. In its three dimensional form, the FBP takes the 
following form: 

™ e -(N+jp)r 

V(r) = -Y,°i " • (28) 

3=0 

In order to suppress the singular point at r — 0, we take the following conditions: 
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V(0) = constant, 
dV(0) = d 2 V(0) = = d n V(Q) = Q 
dr dr 2 dr n ' ^ ' 

So we can determine a.j by solve the following equations, inferred from the flat bottom condition Eq.(29), 



3=0 



^ aj (N + jp) = V(0), 



3=0 



£>(iV + ip) 2 = 0, 



3=0 



J2 aj (N + j P r = 0. (30) 

3=0 



As in previous literature [Q, n is set to b 9. 

In the rainbow approximation, the SD equation takes the following form: 



S- 1 (p) = l-p-m+-^- J J^L lfl S(k) lv G^(k-p), (31) 



where 



S-'ip) = A(p 2 )j ■ p - B(p 2 ) ee A(p 2 )[ 7 ■ p - m(p 2 )], (32) 
= -0T - ^)G(fc 2 ), (33) 

and m stands for an explicit quark mass-breaking term. With the explicit small mass term, we can preclude the zero 
solution for B(p) and in fact there indeed exist a bare current quark mass. Here we take Landau gauge. This dressing 
comprises the notation of constituent quarks by providing a mass m(p 2 ) = B(p 2 )/A(p 2 ), which is corresponding to 
dynamical symmetry breaking. Because the form of the gluon propagator g 2 G(jp) in the infrared region is unknown 
(as we have discussed in section 2), one often uses model forms as input in the previous studies of the rainbow SD 



equation |29 1 



In this article, we assume that a Wick rotation to Euclidean variables is allowed, and perform a rotation analyti- 
cally continuing p and k into the Euclidean region where them can be denoted by p and k, respectively. Alternatively, 
one can derive the SD equation from the Euclidean path-integral formulation of the theory, thus avoiding possible 
difficulties in performing the Wick rotation . As far as only numerical results are concerned, the two procedures 
are equal. 

The Euclidean SD equation can be projected into two coupled integral equations for A(p 2 ) and B(p 2 ). For sim- 
plicity, here ignore the bar on p and k. Numerical values for A(p 2 ), B{p 2 ) and m(p 2 ) are shown in Fig.[l] . 



VI. CONCLUSION AND DISCUSSION 

In this article, we calculate the dressed quark propagator with the FBP in the framework of the rain-bow SD 
equation. The dressed quark propagator exhibits a dynamical symmetry breaking phenomenon and gives a constituent 
quark mass about 490 MeV (m(m 2 )), which is larger then the value of commonly used constituent quark mass 350 



8 



MeV in the chiral quark model. Then based on dressed quark propagator, we obtain the quark condensate (0\qq\0), the 
7r decay constant. The calculated f n and (0\qq\0) are compatible with experimental and theoretical works respectively. 
Though the value of f w , we can fit the parameters: N = 1.0 x 200MeV, V(0) = -27.0 x 200MeV, p = 5.0 x 200MeV, 
rn u = rrid = 8MeV and the large momentum cut-off L = 630 x 200MeV in solving the SD equation. 

The 7r decay constant /„• is an important parameter in describing the interplay between dynamical symmetry 
breaking and confinement, while the quark condensate serves as an order parameter for dynamical chiral symmetry 
breaking. In the chiral limit, the n plays a double role both as Nambu-Goldstone boson and qq bound state, which 
lead to the result that the scalar part of the inverse dressed quark propagator coincides with the approximated 7r's 
Bethe-Salpeter function. 

Till now, all theoretical studies based on QCD are suffered from shortcomings in one or another way jy], so they 
can not give satisfactory infrared behaviour for gluon propagator which may lead to a linear or oscillator potential 
implements confinement. In this article, we take an alternative description by studying the analytic structure of 
the quark propagator. If the dressed quark propagator is free singularities on the real timelike p 2 axial, although 
the attraction between quarks may diminish at a characteristic distance, but the vacuum repulsive force due to the 
absence of mass poles compresses the quarks together. The traditional confinement based on the arising potential 
between quarks as the separation becomes long can be realized as combination of both the attractive and repulsive 
forces. In our calculation, we take a Fourier transform with respect Euclidean time T, which is often used in lattice 
calculation to extrapolate physical masses, and find that at large time T, there does not exist a physical quark mass, 
which may preclude the existence of physical free quark spectrum. 

In calculation, the coupled integral equations for the quark propagator functions A(p 2 ) and B(p 2 ) are solved 
numerically by simultaneous iterations. The iterations converge rapidly to a unique stable solution of propagator 
functions and independent the initial guesses. The propagator functions A(p 2 ) and B(p 2 ) are shown in Fig.fl] , at 
small p 2 , A(p 2 ) differs from the value 1 appreciably, while it tends to 1 for large p 2 . We find that at small p 2 , m(p 2 ) 
is greatly re-normalized, while at large p 2 , it takes asymptotic behaviour. For u and d quark, m(0) = 632MeV, the 
connection of m(p) to constituent masses is somewhat less direct and is precise only for heavy quarks. For heavy 
quarks, m consUtuent (p) = m(p = 2m constituent (p)) , for light quarks , it only makes a crude estimation At about 
p = IGeV, the mass function grows rapidly as the momentum decreases, that is an indication of dynamical symmetry 
breaking. There is a popularly used constituent quark definition, m(p 2 ) 2 = p 2 ]29| . Here p 2 is Euclidean momentum 
square and m{m 2 ) = 490MeV. Our result is larger than the usually used constituent quark mass m = 350MeV, 
however, it gives a good description of dynamical symmetry breaking. 

The phcnomcnological FBP can be extended to other system, such as qq system, qQ system, vector meson and 
qqq system or by applying the results to the calculation of quantities and processes requiring detailed knowledge of 
the quark propagator . 
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